An apparatus and a method for generating data representative of a pixel beam

ABSTRACT

There are several types of plenoptic devices having their proprietary file format. At present there is no standard supporting the acquisition and transmission of multidimensional information for an exhaustive over-view of the different parameters upon which a light-field depends. As such acquired light-field data for different cameras have a diversity of formats. The notion of pixel beam, which represents a volume occupied by a set of rays of light in an object space of an optical system of a camera is thus introduced. The method according to the invention enables to provide data representative of a collection of pixel beams describing a first optical system that is agnostic since these data are obtained by imaging the collection of pixel beams through a second optical system. Such data representative of a collection of pixel beams enable the generation of parametrized output images from which post-processing.

TECHNICAL FIELD

The present invention relates to generation of data representing a light field.

BACKGROUND

The acquisition of four-dimensional or 4D light-field data, which can be viewed as a sampling of a 4D light field, i.e. the recording of light rays, is explained in the article “Understanding camera trade-offs through a Bayesian analysis of light field projections” by Anat Levin and al., published in the conference proceedings of ECCV 2008 is an hectic research subject.

Compared to classical two-dimensional or 2D images obtained from a camera, 4D light-field data enable a user to have access to more post-processing features that enhance the rendering of images and the interactivity with the user. For example, with 4D light-field data, it is possible to perform refocusing of images with freely selected distances of focalization meaning that the position of a focal plane can be specified/selected a posteriori, as well as changing slightly the point of view in the scene of an image. In order to acquire 4D light-field data, several techniques can be used. For example, a plenoptic camera is able to acquire 4D light-field data. A plenoptic camera comprises a main lens, a microlens array comprising a plurality of micro-lenses arranged in a two-dimensional array and an image sensor. Another way to acquire 4D light-field data is to use a camera array which comprises a lens array and an image sensor.

In the example of the plenoptic camera, the main lens receives light from an object in an object field of the main lens and passes the light through an image field of the main lens.

At last, another way of acquiring a 4D light field is to use a conventional camera that is configured to capture a sequence of 2D images of a same scene at different focal planes. For example, the technique described in the document “Light ray field capture using focal plane sweeping and its optical reconstruction using 3D displays” by J.-H. Park et al., published in OPTICS EXPRESS, Vol. 22, No. 21, in October 2014, may be used to achieve the acquisition of 4D light field data by means of a conventional camera.

There are several ways to represent 4D light-field data. Indeed, in the Chapter 3.3 of the Ph. D dissertation thesis entitled “Digital Light Field Photography” by Ren Ng, published in July 2006, three different ways to represent 4D light-field data are described. Firstly, 4D light-field data can be represented, when recorded by a plenoptic camera by a collection of micro-lens images. 4D light-field data in this representation are named raw images or raw 4D light-field data. Secondly, 4D light-field data can be represented, either when recorded by a plenoptic camera or by a camera array, by a set of sub-aperture images. A sub-aperture image corresponds to a captured image of a scene from a point of view, the point of view being slightly different between two sub-aperture images. These sub-aperture images give information about the parallax and depth of the imaged scene. Thirdly, 4D light-field data can be represented by a set of epipolar images see for example the article entitled: “Generating EPI Representation of a 4D Light Fields with a Single Lens Focused Plenoptic Camera”, by S. Wanner and al., published in the conference proceedings of ISVC 2011.

Light-field acquisition devices are extremely heterogeneous. Light-field cameras are of different types for example plenoptic or camera arrays. Within each type there are many differences such as different optical arrangements, or micro-lenses of different focal lengths. Each camera has its own proprietary file format. At present here is no standard supporting the acquisition and transmission of multi-dimensional information for an exhaustive over-view of the different parameters upon which a light-field depends. As such acquired light-field data for different cameras have a diversity of formats.

The present invention has been devised with the foregoing in mind.

SUMMARY OF INVENTION

According to a first aspect of the invention there is provided a computer implemented method for generating data representative of a volume occupied by a set of rays of light passing through a pupil of an optical system, said volume occupied by said set of rays of light being called a pixel beam, comprising:

-   -   acquiring data representative of said pixel beam in an object         space of a first optical system, said pixel beam being occupied         by a set of rays of light passing through a pupil of said first         optical system and an image conjugate, in said object space of         the first optical system, of at least one pixel of a sensor         associated to said first optical system;     -   generating data representative of said pixel beam by computing         an image conjugate of said pixel beam through a second optical         system. Such a method enables to provide data representative of         a collection of pixel beams describing a first optical system         that is agnostic, i.e. which is independent of the camera         embedding the optical system, since these data are obtained by         imaging the collection of pixel beams through a second optical         system. Such data representative of a collection of pixel beams         enable the generation of parametrized output images from which         post-processing can be performed such as creating focal stacks,         refocusing, perspective changes, mixing of light-field contents         acquired from different acquisition systems, etc.

One of the goals of computing an image of a collection of pixel beams describing a first optical system through a second optical system is to provide a representation of the object space of several optical systems embedded in various cameras in a dual and compact space which corresponds to a focal volume of the optical system. Such a focal volume can be easily projected on an image plane for example during a refocusing operation.

An advantage of imaging a collection of pixel beams describing a first optical system through a second optical system is that the imaged pixel beams are located in a single focal volume and that collections of pixel beams representing various types of optical acquisition devices can be positioned in a same focal volume using the same optical system which make processing of the data representative of the different collections of pixel beams easier.

According to another embodiment of the invention, the method further comprises:

-   -   computing an intersection, in the object space of the first         optical system, of said pixel beam with a plane located at a         distance z from a centre of the pupil,     -   computing an image conjugate of said intersection through said         second optical system to provide data representative of the         pixel beam.

Such a method enables to perform synthetic aperture refocusing. Synthetic aperture imaging consists in rearranging images acquired by an optical system embedded in a camera in order to simulate another camera for refocusing purposes. The pixel beams representing light-field data acquired by a first camera are thus rearranged so that objects belonging to a given plane are in focus. The plane is defined by its affine equation in a reference coordinate system.

According to another embodiment of the invention, the first optical system is an actual camera and the second optical system is a virtual camera.

In such an embodiment of the invention, the knowledge of the properties of the entrance pupils of the first optical system is required and of the properties of the exit pupils of the second optical system in order to compute the conjugate of the pixel beam in the image space of the second optical system.

As long as the second optical system, i.e. the virtual camera, images precisely its focus plane, the precise knowledge of the properties of its exit pupils is not required since the image of the pixel beam through the second optical system is located on the virtual sensor.

According to another embodiment of the invention, the method the first optical system is a virtual camera and the second optical system is an actual camera.

In such an embodiment of the invention, the knowledge of the properties of the entrance pupils of the first optical system is required and of the properties of the exit pupils of the second optical system, i.e. the actual camera, in order to compute the conjugate of the pixel beam in the image space of the second optical system.

As long as the first optical system, i.e. the virtual camera, images precisely its focus plane, the precise knowledge of the properties of its exit pupils is not required since the pixel beam are the object conjugates of the pixels of a virtual sensor associated with the first acquisition system.

Another object of the invention concerns a device for rendering an image from light field data using obtained in accordance with the method according to different embodiments of the invention.

Another object of the invention concerns a device for generating data representative of a volume occupied by a set of rays of light passing through a pupil of an optical system, said volume occupied by said set of rays of light being called a pixel beam, the device comprising a light field data generation module configured to:

-   -   acquire data representative of said pixel beam in an object         space of a first optical system, said pixel beam being occupied         by a set of rays of light passing through a pupil of said first         optical system and an image conjugate, in said object space of         the first optical system, of at least one pixel of a sensor         associated to said first optical system;     -   generate data representative of said pixel beam by computing an         image conjugate of said pixel beam through a second optical         system.

Another object of the invention concerns a light field imaging device comprising:

-   -   an array of micro lenses arranged in a regular lattice         structure;     -   a photosensor configured to capture light projected on the         photosensor from the array of micro lenses, the photosensor         comprising sets of pixels, each set of pixels being optically         associated with a respective micro lens of the array of micro         lenses; and     -   a device for generating metadata in accordance with claim 8.

Another object of the invention concerns a data package for data representative of a volume in an object space of a first optical system occupied by a set of rays of light passing through a pupil of said first optical system and a conjugate, in said object space of the first optical system, of at least one pixel of a sensor associated to said first optical system, said volume occupied by said set of rays of light being called a pixel beam said data representative of a pixel beam being an image conjugate of said pixel beam through a second optical system.

Some processes implemented by elements of the invention may be computer implemented. Accordingly, such elements may take the form of an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining software and hardware aspects that may all generally be referred to herein as a “circuit”, “module” or “system’. Furthermore, such elements may take the form of a computer program product embodied in any tangible medium of expression having computer usable program code embodied in the medium.

Since elements of the present invention can be implemented in software, the present invention can be embodied as computer readable code for provision to a programmable apparatus on any suitable carrier medium. A tangible carrier medium may comprise a storage medium such as a floppy disk, a CD-ROM, a hard disk drive, a magnetic tape device or a solid state memory device and the like. A transient carrier medium may include a signal such as an electrical signal, an electronic signal, an optical signal, an acoustic signal, a magnetic signal or an electromagnetic signal, e.g. a microwave or RF signal.

BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the invention will now be described, by way of example only, and with reference to the following drawings in which:

FIG. 1 shows a volume occupied by a set of rays of light passing through a pupil of an optical system 11 of a camera,

FIG. 2 shows a hyperboloid of one sheet representing a pixel beam according to an embodiment of the invention,

FIG. 3 shows a hyperboloid of one sheet and its asymptotic cones,

FIG. 4 shows a pixel beam represented by two coaxial, partially overlapping cones according to an embodiment of the invention,

FIG. 5 is a schematic block diagram illustrating an example of an apparatus for generating data representative of pixel beams according to an embodiment of the invention,

FIG. 6 is a flow chart for explaining a process for encoding an image captured by the optical acquisition system according to an embodiment of the disclosure,

DETAILED DESCRIPTION

As will be appreciated by one skilled in the art, aspects of the present principles can be embodied as a system, method or computer readable medium. Accordingly, aspects of the present principles can take the form of an entirely hardware embodiment, an entirely software embodiment, (including firmware, resident software, micro-code, and so forth) or an embodiment combining software and hardware aspects that can all generally be referred to herein as a “circuit”, “module”, or “system”. Furthermore, aspects of the present principles can take the form of a computer readable storage medium. Any combination of one or more computer readable storage medium (a) may be utilized.

For any optical acquisition system, may it be plenoptic or not, in addition to raw images or epipolar images representing 4D light-field data captured by the optical acquisition system, it is interesting to obtain information related to a correspondence between pixels of a sensor of said optical acquisition system and an object space of said optical acquisition system. Knowing which portion of the object space of an optical acquisition system a pixel belonging to the sensor of said optical acquisition system is sensing enables the improvement of signal processing operations such as de-multiplexing, de-mosaicking, refocusing, etc., and the mixing of images captured by different optical systems with different characteristics. Furthermore, information related to the correspondence between the pixels of the sensor of the optical acquisition system and the object space of said optical acquisition system are independent of the optical acquisition system.

The present disclosure introduces the notion of pixel beam 10, shown on FIG. 1, which represents a volume occupied by a set of rays of light passing through a pupil of an optical system 11 of a camera (not shown on FIG. 1) and a conjugate of a pixel of a sensor of the camera in an object space of the optical system in a direction normal to a surface of the pupil

The set of rays of light is sensed by a pixel 12 of a sensor 13 of the camera through a pupil 14 of said optical system 11. The optical system 11 may be a combination of lenses fit for photo or video cameras. A pupil of an optical system is defined as the image of an aperture stop as seen through said optical system, i.e. the lenses of the optical acquisition system, which precedes said aperture stop. An aperture stop is an opening which limits the amount of light which passes through the optical system of the optical acquisition system. For example, an adjustable blade diaphragm located inside a camera lens is the aperture stop for the lens. The amount of light admitted through the diaphragm is controlled by the diameter of the diaphragm opening which may adapted depending of the amount of light a user of the camera wishes to admit. For example, making the aperture smaller reduces the amount of light admitted through the diaphragm, and, simultaneously, increases the depth of focus. The apparent size of a stop may be larger or smaller than its physical size because of the refractive action of a portion of the lens. Formally, a pupil is the image of the aperture stop through all lenses of the optical acquisition system located between the physical stop and the observation space.

A pixel beam 10 is defined as a pencil of rays of light that reach a given pixel 12 when propagating through the optical system 11 via an entrance pupil 14. As light travels on straight lines in free space, the shape of such a pixel beam 10 can be defined by two sections, one being the conjugate 15 of the pixel 12, and the other being the entrance pupil 14. The pixel 12 is defined by its non-null surface and its sensitivity map.

In a first embodiment of the invention, a pixel beam 30 may be represented by an hyperboloid of one sheet, as shown on FIG. 2, supported by two elements: the pupil 24 and the conjugate 25 of the pixel 12 in the object space.

A hyperboloid of one sheet is a ruled surface that can support the notion of pencil of rays of light and is compatible with the notion of “étendue” of physical light beams, notion linked to the preservation of energy across sections of the physical light beams.

As represented on FIG. 3, a hyperboloid of one sheet 30 is mostly identical to its asymptotic cones 31, 32, except in the fundamental region of its smallest section, called the waist 35, which corresponds to the conjugate 15 in the object space. For plenoptic systems, such as light-field cameras, this is the region where space sampling by multiple path rays is performed. Sampling space with unique cones degenerating to a point in this region is not adequate, as pixel 12 sensitivity is significant on some tens of square microns on its surface and cannot be represented by a mathematical point with infinitely small surface as would be a cone tip.

In an embodiment of the invention, each pixel beam 10, 20, 30, is defined by four independent parameters: z_(P), θ_(x), θ_(y), a defining the position and size of the pixel conjugate 15, 35, in front of the pupil 14, 24 and by six pupilar parameters x_(O), y_(O), z_(O), θ_(x0), θ_(y0), r which define the position, orientation and radius of the pupil 14, 24. These six pupilar parameters are common to the collection of pixel beams, when represented by a hyperboloid of one sheet, sharing a same pupil 14, 24. Indeed, a pixel beam represents the volume occupied by a set of rays of light in the object space of the optical system 11 sensed by the pixel 12 through the pupil 14, i.e. to a given couple pixel 12/pupil 14, 24 corresponds a unique pixel beam 10, 20, 30, but a plurality of distinct pixel beams can be supported by a same pupil 14, 24.

An origin O of a coordinate system (x, y, z) in which the parameters of the hyperboloid of one sheet representing the pixel beam 10, 20, 30 are defined corresponds to the centre of the pupil 14 as shown on FIG. 1, where the z axis defines a direction normal to the surface of the pupil 14, 24.

The parameters θ_(x), θ_(y), define chief ray directions relative to the entrance of the pupil 14 centre. They depend on the pixel 12 position on the sensor 13 and on the optical elements of the optical system 11. More precisely, the parameters θ_(x), θ_(y) represent shear angles defining a direction of the conjugate 15 of the pixel 12 from the centre of the pupil 14.

The parameter z_(P) represents a distance of the waist 35 of the pixel beam 10, 20, 30, or the conjugate 15 of the pixel 12, along the z axis.

The parameter a represents the radius of the waist 35 of the pixel beam 10, 20, 30.

For optical systems 11 where optical distortions and field curvatures may be modelled, the parameters z_(P) and a can depend on the parameters θ_(x) and θ_(y) via parametric functions.

The four independent parameters are related to the pixel 12 and its conjugate 15.

The six complementary pupilar parameters defining a pixel beam 10, 20, 30 are:

-   -   r which represents the pupil 14, 24 radius,     -   x_(O), y_(O), z_(O) which represent the coordinates of the pupil         14, 24 centre in the (x, y, z) coordinate system, and     -   θ_(x0), θ_(y0) which represent the orientation of the pupil 14,         24 in the reference (x, y, z) coordinate system.

These six pupilar parameters are related to the pupil 14, 24. Another parameter c is defined. Such a parameter c is dependent on the parameters z_(P) and a related to the pixel 12 and its conjugate 15 and on the parameters r related to the pupil 14, 24. The parameter c defines the angular aperture α of the pixel beam 10, 20, 30 and is given by the formula:

${\tan (\alpha)} = {\frac{a}{c}.}$

Thus the expression of the parameter c is given by the following equation:

$\begin{matrix} {c^{2} = \frac{a^{2}z_{P}^{2}}{r^{2} - a^{2}}} & (1) \end{matrix}$

The coordinates (x, y, z), in the object space, of points belonging to the surface delimiting the pixel beam 10, 20, 30 are function of the above defined sets of parameters related to the pupil 14, and to the conjugate 15 of the pixel. Thus, equation (2) enabling the generation of the hyperboloid of one sheet representing the pixel beam 10, 20, 30 is:

$\begin{matrix} {{\frac{\left( {x - {z \cdot {\tan \left( \theta_{x} \right)}}} \right)^{2}}{a^{2}} + \frac{\left( {y - {z \cdot {\tan \left( \theta_{y} \right)}}} \right)^{2}}{a^{2}} - \frac{\left( {z - z_{P}} \right)^{2}}{c^{2}}} = 1} & (2) \end{matrix}$

A parametric equation (3) of the same hyperboloid representing the pixel beam 10, 20, 30 is:

$\begin{matrix} \left\{ \begin{matrix} {x = {{a{\sqrt{1 + \frac{\left( {z - z_{P}} \right)^{2}}{c^{2}}} \cdot {\cos (v)}}} + {z \cdot {\tan \left( \theta_{x} \right)}}}} \\ {y = {{a{\sqrt{1 + \frac{\left( {z - z_{P}} \right)^{2}}{c^{2}}} \cdot {\sin (v)}}} + {z \cdot {\tan \left( \theta_{y} \right)}}}} \end{matrix} \right. & (3) \end{matrix}$

wherein ν is an angle in the (x, y) plane enabling the generation of the pixel beam 10, 20, 30 from a generating hyperbola, ν varies in [0, 2π] interval, and z∈[0, ∞] is the coordinate along the z axis which defines a direction normal to the surface of the pupil 14, 24. Equations (2) and (3) are written on the assumption that the section of the pixel 12 and its conjugate 15 are circular and that the section of the pupil 14, 24 is circular as well.

Information related to a correspondence between pixels of a sensor of said optical acquisition system and an object space of said optical acquisition system may take the form of either a set of parameters comprising the four independent parameters: z_(P), θ_(x), θ_(y), a defining the position and size of the pixel conjugate 15, 35, in front of the pupil 14, 24 and the six pupilar parameters x_(O), y_(O), z_(O), θ_(x0), θ_(y0), r which define the position, orientation and radius of the pupil 14, 24 when the pixel beam is to be represented by its parametric equation.

Thus, this set of parameters is provided in addition to raw images or epipolar images representing 4D light-field data captured by the optical acquisition system in order to be used while processing the 4D light-field data.

In a second embodiment of the invention, a pixel beam 40 may be represented by two coaxial, partially overlapping cones a front cone 41 _(F) and a rear cone 41 _(R) as shown on FIG. 4, supported by two elements: the pupil 44 and the conjugate 45 of the pixel 42 in the object space, i.e. the surface in the object space that is imaged on the pixel.

The front cone 41 _(F) is the image of a convex frustum defined by the pixel 42 and the pupil 44. The apex of the convex frustum lies beyond the sensor of the optical acquisition system. By construction, the front cone 41 _(F) is converging in the object space of the optical acquisition system and the apex of the front cone 41 _(F) lies between the conjugate of the pixel 45, or the waist of the pixel beam 40, and the pupil dd. The front cone 41 _(F) derives from the solid angle subtended by the pupil 44 at the pixel 42.

The rear cone 41 _(R) is the image of a cone defined by the pixel 42 and the pupil 44, the apex of which lies between the pupil 44 and the sensor of the optical acquisition system. By construction, the apex of the rear cone 41 _(R) is located beyond the waist 45 of the pupil 40. The rear cone 41 _(R) does not necessarily converge in the object space of the optical acquisition system, in some cases, it may degenerate into a cylinder or a diverging cone. In the latter case, the apex of the diverging cone lies in the image space of the optical acquisition system, i.e. before the entrance of the pupil 44.

The front cone 41 _(F) and the rear cone 41 _(R) share the same revolution axis, which is a line joining the centre of the pupil 44 and the centre of the waist 45.

Cones are ruled surfaces that can support the notion of pencil of rays of light and when combining two cones is compatible with the notion of “étendue” of physical light beams, notion linked to the preservation of energy across sections of the physical light beams. Intersections of cones with planes are conic curves, as for hyperboloids, which can be characterized by a plurality of coefficients. Considering its apex, a cone may be represented by three angular parameters: a polar angle measured from the revolution axis of the cone, up to the apex angle and the direction of the revolution axis given by two angles.

Let xyz be the coordinate system of the optical acquisition system, z denoting the optical axis of the optical acquisition system with z>0 in the object space of the optical acquisition system and the centre of the pupil 44 being the origin of said coordinate system.

The optics of the optical acquisition system images the object space of the optical acquisition system from the range z∈[2f; +∞] into the image space of the optical acquisition system z∈[−2f; −f], where f is the focal length of the optics of the optical acquisition system. The location of the pupil 44 and the waist 45 of the pixel beam 40 are known in the coordinate system xyz of the optical acquisition system from the calibration of the optical acquisition system. The pupil 44 and the waist 45 are assumed to be parallel and are both normal to the z axis.

Let us call z′ the chief ray of the pixel beam 40. The chief ray is the line joining the centre of the pupil 44 and the centre of the waist 45 of the pixel beam 40. The chief ray is also the revolution axis and the axis of symmetry of the pixel beam 40. Thus, in the coordinate system xyz′, the pixel beam 40 is a solid of revolution.

Both the apices of the front cone 41 _(F) and the rear cone 41 _(R) are located on the chief ray z′ of the pixel beam 40. Under the thin lens approximation, the coordinates of these two apices are computed in the coordinate system xyz of the optical acquisition system as follow, under the assumption that the sensor of the optical acquisition system is not located the rear focal plane:

$\begin{matrix} \left\{ \begin{matrix} {\frac{z_{w} - z_{front}}{W} = \frac{z_{front} - z_{P}}{P}} \\ {\frac{z_{rear} - z_{w}}{W} = \frac{z_{rear} - z_{P}}{P}} \end{matrix} \right. & (4) \end{matrix}$

i.e.:

$\begin{matrix} \left\{ \begin{matrix} {z_{front} = \frac{{P \cdot z_{W}} + {W \cdot z_{P}}}{P + W}} \\ {z_{rear} = \frac{{P \cdot z_{W}} - {W \cdot z_{P}}}{P - W}} \end{matrix} \right. & (5) \end{matrix}$

where P, z_(P), W and z_(w) respectively denote the diameter of the pupil 44 with P>0, its z-coordinate, the diameter of the pixel's conjugate 45 with 0<W<+∞, and its z-coordinate 0<z_(w)<+∞.

The z-coordinate z_(rear) of the apex of the rear cone 41 _(R) may be positive, when the rear cone 41 _(R) is a converging cone, negative, when the rear cone 41 _(R) is a diverging cone. It may also be infinite if the pupil 44 and the pixel's conjugate 45 of the pixel beam are of the same size.

If the sensor of the optical acquisition system is located on the rear focal plane, then W=+∞ and z_(w)=+∞. As their ratio is a constant:

$\begin{matrix} {\frac{1}{z_{front} - z_{P}} = {\frac{p}{P \cdot f} = \frac{1}{z_{P} - z_{rear}}}} & (6) \end{matrix}$

where p and f respectively represent the diameter of the pixel 42 with p>0 and the focal length of the optics of the optical acquisition system with f>0 assuming the optics of the optical acquisition system is a converging lens.

-   -   The apex angles are given by:

$\begin{matrix} \left\{ \begin{matrix} {{\tan \; \theta_{front}} = \frac{W/2}{{z_{front} - z_{w}}}} \\ {{\tan \; \theta_{rear}} = \frac{W/2}{{z_{rear} - z_{w}}}} \end{matrix} \right. & (7) \end{matrix}$

-   -   Considering the apex of each cones, which union represents a         pixel beam 40, rays can be defined with two angular parameters:         the polar angle measure from the revolution axis of the pixel         beam, up to the apex angle, and an azimuth in [0,2π[.

Those information related to pixel beams are metadata associated to a given optical acquisition system. They may be provided as a data file stored for example on a CD-ROM or a flash drive supplied with the optical acquisition system. The data file containing the additional information related to pixel beams may also be downloaded from a server belonging to the manufacturer of the optical acquisition system. In an embodiment of the invention, these additional information related to pixel beams may also be embedded in a header of the images captured by the optical acquisition system.

The knowledge of these information related to pixel beams enables the processing of images captured by any optical acquisition system independently of the proprietary file format and of the features of the optical acquisition system used to capture the images to be processed.

The knowledge of information related to pixel beams enables the processing of images captured by any optical acquisition system independently of the proprietary file format and of the features of the optical acquisition system used to capture the images to be processed.

FIG. 5 is a schematic block diagram illustrating an example of an apparatus for generating data representative of pixel beams representing the object space of a first optical system by imaging said pixel beams through a second optical system according to an embodiment of the present disclosure.

The apparatus 500 comprises a processor 501, a storage unit 502, an input device 503, a display device 504, and an interface unit 505 which are connected by a bus 506. Of course, constituent elements of the computer apparatus 500 may be connected by a connection other than a bus connection.

The processor 501 controls operations of the apparatus 500. The storage unit 502 stores at least one program capable of generating data representative of pixel beams representing the object space of a first optical system when these pixel beams are imaged through a second optical system to be executed by the processor 501, and various data, including parameters related to a position of the pixel 12 on the sensor 13 or parameters related to the first optical system 11 of the optical acquisition system and a second optical system, parameters used by computations performed by the processor 501, intermediate data of computations performed by the processor 501, and so on. The processor 501 may be formed by any known and suitable hardware, or software, or a combination of hardware and software. For example, the processor 501 may be formed by dedicated hardware such as a processing circuit, or by a programmable processing unit such as a CPU (Central Processing Unit) that executes a program stored in a memory thereof.

The storage unit 502 may be formed by any suitable storage or means capable of storing the program, data, or the like in a computer-readable manner. Examples of the storage unit 502 include non-transitory computer-readable storage media such as semiconductor memory devices, and magnetic, optical, or magneto-optical recording media loaded into a read and write unit. The program causes the processor 501 to perform a process for computing data representative of the pixel beams of a collection of pixel beams representative of the object space of a first optical system from an image conjugate of said pixel beam through a second optical system beam according to an embodiment of the present disclosure as described hereinafter with reference to FIG. 6.

The input device 503 may be formed by a keyboard, a pointing device such as a mouse, or the like for use by the user to input commands, to make user's selections of parameters used for generating a parametric representation of a volume occupied by a set of rays of light in an object space of an optical system. The output device 604 may be formed by a display device to display, for example, a Graphical User Interface (GUI), images generated according to an embodiment of the present disclosure. The input device 503 and the output device 504 may be formed integrally by a touchscreen panel, for example.

The interface unit 505 provides an interface between the apparatus 500 and an external apparatus. The interface unit 505 may be communicable with the external apparatus via cable or wireless communication. In an embodiment, the external apparatus may be an optical acquisition system such as an actual camera.

FIG. 6 is a flow chart for explaining a process for encoding an image captured by the optical acquisition system according to an embodiment of the invention.

The processor 501 of the apparatus 500 executes the program capable of computing data representative of at least one pixel beam 70 of a collection of pixel beams representative of the object space of a first optical system (not shown on the figure) from an image conjugate 72 of said pixel beam 70 through a second optical system 71 as represented on FIG. 7. In the following example, the pixel beams 70, 72 are represented by a hyperboloid of one sheet 20.

Thus, the processor 501 computes an image conjugate 72 of the pixel beam 70 through the optical system 71, which is for example a thin lens, to provide data representative of the pixel beam 70 in the focal volume of the optical system 71.

The image conjugate 72 of the pixel beam 70 through the optical system 71 may be computed as explained hereafter.

As described in reference to FIGS. 1 to 3, a pixel beam linked to a pupil of an optical system is defined by six extrinsic parameters:

-   -   r: the entrance pupil radius of the first optical system,     -   x_(O), y_(O), z_(O): the coordinates of the entrance pupil         centre of the first optical system in a reference coordinate         system, and     -   θ_(x0), θ_(y0): an orientation of the entrance pupil of the         first optical system in the reference coordinate system.

In the reference coordinate system, and keeping θ_(x0)=θ_(y0)=0 for simplicity, a parametric equation of the hyperboloid representing the pixel beam is given by:

$\begin{matrix} {{\frac{\left( {x - x_{O} - {\left( {z - z_{0}} \right) \cdot {\tan \left( \theta_{x} \right)}}} \right)^{2}}{a^{2}} + \frac{\left( {y - y_{O} - {\left( {z - z_{0}} \right) \cdot {\tan \left( \theta_{y} \right)}}} \right)^{2}}{a^{2}} - \frac{\left( {z - z_{O} - z_{P} - z_{O}} \right)^{2}}{c^{2}}} = 1} & (8) \end{matrix}$

or equivalently:

$\begin{matrix} \left\{ \begin{matrix} {x = {x_{O} + {a{\sqrt{1 + \frac{\left( {z - z_{P}} \right)^{2}}{c^{2}}} \cdot {\cos (v)}}} + {\left( {z - z_{0}} \right) \cdot {\tan \left( \theta_{x} \right)}}}} \\ {y = {y_{O} + {a{\sqrt{1 + \frac{\left( {z - z_{P}} \right)^{2}}{c^{2}}} \cdot {\sin (v)}}} + {\left( {z - z_{0}} \right) \cdot {\tan \left( \theta_{y} \right)}}}} \end{matrix} \right. & (9) \end{matrix}$

In an embodiment of the invention, a value representing a colour, such as a RGB value, acquired by a pixel of the sensor of the first optical acquisition system is associated to the corresponding pixel beam 70.

The parameters defining a plane surface representing the second optical system 71 are either defined when the second optical system 71 is a virtual optical system or retrieved when the second optical system is an actual optical system:

-   -   coordinates of the optical centre of the second optical system         71 are x_(C), y_(C), z_(C)     -   coordinates of a unit vector normal to the surface of the         optical system 71: θ_(xC), θ_(yC)     -   a focal f of the optical system 71.     -   In order to simplify explanations, the optical centre of the         second optical system 71 is taken in the same plane as the pupil         centre of the first optical system: z_(C)=z_(O) and the optical         system 71 is supposed to be parallel to the pupil, i.e.         θ_(xC)=θ_(yC)=0.

The coordinates of a point M relative to the optical system 71 center C are denoted: x, y, z. The coordinates of the image M′ of point M through the optical system 71, are denoted: x′, y′, z′. Each of these coordinates are algebraic values.

-   -   The coordinates of point M′ are computed using the Newton's form         of “the Lens Maker formula” which is:

(|z|−f)(|z′|−f)=f ²  (10)

Since point M is located before the optical system 71, i.e. z<0 and point M′ is located after the optical system 71, i.e. z′>0, then:

(−z−f)(z′−f)=f ²  (11)

which gives

$\begin{matrix} {z^{\prime} = {{f - \frac{f^{2}}{z + f}} = \frac{f \cdot z}{z + f}}} & (12) \end{matrix}$

When the optical system 71 is a thin lens, as the ray of light passing through its optical centre is not deviated magnification is given by:

$m = {\frac{x^{\prime}}{x} = {\frac{y^{\prime}}{y} = \frac{z^{\prime}}{z}}}$

which gives:

$\begin{matrix} \left\{ \begin{matrix} {z^{\prime} = \frac{f \cdot z}{f + z}} \\ {x^{\prime} = {x\frac{z^{\prime}}{z}}} \\ {y^{\prime} = {y\frac{z^{\prime}}{z}}} \end{matrix} \right. & (13) \\ {{and}\mspace{14mu} {conversely}\text{:}} & \; \\ \left\{ \begin{matrix} {z = \frac{f \cdot z^{\prime}}{f - z^{\prime}}} \\ {x = {x^{\prime}\frac{z}{z^{\prime}}}} \\ {y = {y^{\prime}\frac{z}{z^{\prime}}}} \end{matrix} \right. & (14) \end{matrix}$

When M is a point belonging to the surface of a pixel beam 70 the coordinates of point M′ situated on the surface of the imaged pixel beam are:

$\begin{matrix} \left\{ \begin{matrix} {z^{\prime} = \frac{f \cdot z}{f + z}} \\ {x^{\prime} = {x\frac{z^{\prime}}{z}}} \\ {y^{\prime} = {y\frac{z^{\prime}}{z}}} \end{matrix} \right. & (15) \end{matrix}$

As a paraxial lens transforms rays of light (straight lines) in object space into rays of light in image space, and as pixel beams 70 are represented by a ruled surface, the imaged pixel beams 72 are also represented by a ruled surface, as the image of an hyperboloid of one sheet in object space is an hyperboloid of one sheet in image space.

Let's determine the parameters of the hyperboloid representing the imaged pixel beam 72 from the parameters of the hyperboloid representing the pixel beam 70 and from the parameters of the optical system 71.

Considering the following equation:

$\begin{matrix} \left\{ \begin{matrix} {x = {x_{O} + {a{\sqrt{1 + \frac{\left( {z - z_{P}} \right)^{2}}{c^{2}}} \cdot {\cos (v)}}} + {\left( {z - z_{0}} \right) \cdot {\tan \left( \theta_{x} \right)}}}} \\ {y = {y_{O} + {a{\sqrt{1 + \frac{\left( {z - z_{P}} \right)^{2}}{c^{2}}} \cdot {\sin (v)}}} + {\left( {z - z_{0}} \right) \cdot {\tan \left( \theta_{y} \right)}}}} \end{matrix} \right. & (9) \end{matrix}$

representing the pixel beam 70 where θ_(x)=θ_(y)=0, x₀=y₀=z₀=0 and 84=0 in order to reduce the analysis to a pixel beam on axis, and to the x, z plane.

The parameters to be determined are then z′_(W), a′ which are a function of z_(P), a, f, r, where z′_(W) is the z coordinate of the waist of the imaged pixel beam 72 (W′ is not the conjugate point of P through the imaging process due to beam aperture—a≠0) and a′ is the radius of the waist of the imaged pixel beam 72.

With the assumption that θ_(x)=θ_(y)=0, x₀=y₀=z₀=0 and ν=0 the equation 16 applied to point M simplifies to:

$\begin{matrix} \left\{ \begin{matrix} {x = {a\sqrt{1 + \frac{\left( {z - z_{P}} \right)^{2}}{c^{2}}}}} \\ {y = 0} \end{matrix} \right. & (16) \end{matrix}$

Using the paraxial lens transform of equation 14 for M we get:

$\begin{matrix} {{x^{\prime}\frac{z}{z^{\prime}}} = {a\sqrt{1 + \frac{\left( {\frac{f \cdot z^{\prime}}{f - z^{\prime}} - z_{P}} \right)^{2}}{c^{2}}}}} & (17) \\ {\left. \Leftrightarrow\frac{x^{\prime^{2}}}{a^{2}} \right. = {\left( \frac{f - z^{\prime}}{f} \right)^{2}\left( {1 + \frac{\left( {\frac{f \cdot z^{\prime}}{f - z^{\prime}} - z_{P}} \right)^{2}}{c^{2}}} \right)}} & (18) \\ {\left. \Leftrightarrow\frac{x^{\prime^{2}}}{a^{2}} \right. = {\frac{\left( {f - z^{\prime}} \right)^{2}}{f^{2}} + \frac{\left( {{\left( {f + z_{P}} \right) \cdot z^{\prime}} - {f \cdot z_{P}}} \right)^{2}}{f^{2}c^{2}}}} & (19) \\ {\left. \Leftrightarrow{\frac{f^{2}c^{2}}{a^{2}}x^{\prime^{2}}} \right. = {{c^{2}\left( {f - z^{\prime}} \right)}^{2} + \left( {{\left( {f + z_{P}} \right) \cdot z^{\prime}} - {f \cdot z_{P}}} \right)^{2}}} & (20) \end{matrix}$

which is the equation representing the imaged pixel beam 72 in the focal volume of the optical system 71.

Deriving e equation (20) enables the calculation of the minima and/or maxima of the hyperbolas representing the pixel beams by finding the zero crossing of the derivative.

Equation (20) when developed gives:

$\begin{matrix} {{\frac{f^{2}c^{2}}{a^{2}}x^{\prime^{2}}} = {{\left( {\left( {f + z_{P}} \right)^{2} + c^{2}} \right) \cdot z^{\prime^{2}}} - {2{{f\left( {{z_{P}\left( {f + z_{P}} \right)} + c^{2}} \right)} \cdot z^{\prime}}} + {f^{2}\left( {z_{P}^{2} + c^{2}} \right)}}} & (21) \end{matrix}$

The derivative of equation (21) is given by:

$\begin{matrix} {{\frac{d}{{dz}^{\prime}}\left( {\frac{f^{2}c^{2}}{a^{2}}x^{\prime^{2}}} \right)} = {\frac{d}{{dz}^{\prime}}\left( {{\left( {\left( {f + z_{P}} \right)^{2} + c^{2}} \right)z^{\prime^{2}}} - {2{{f\left( {{z_{P}\left( {f + z_{P}} \right)} + c^{2}} \right)} \cdot z^{\prime}}} + {f^{2}\left( {z_{P}^{2} + c^{2}} \right)}} \right)}} & (22) \\ {\mspace{76mu} {\left. \Leftrightarrow{\frac{d}{{dz}^{\prime}}\left( {\frac{f^{2}c^{2}}{a^{2}}x^{\prime^{2}}} \right)} \right. = {{2\left( {\left( {f + z_{P}} \right)^{2} + c^{2}} \right)z^{\prime}} - {2{f\left( {{z_{P}\left( {f + z_{P}} \right)} + c^{2}} \right)}}}}} & (23) \end{matrix}$

Computing the zero crossing gives of the derivative (23):

$\begin{matrix} {{\frac{d}{{dz}^{\prime}}\left( {\frac{f^{2}c^{2}}{a^{2}}x^{\prime 2}} \right)} = {\left. 0\Leftrightarrow z_{w}^{\prime} \right. = \frac{f\left( {{z_{P}\left( {f + z_{P}} \right)} + c^{2}} \right)}{\left( {f + z_{P}} \right)^{2} + c^{2}}}} & (24) \\ {\left. \Leftrightarrow z_{w}^{\prime} \right. = {f\left( {1 - \frac{f\left( {f + z_{P}} \right)}{\left( {f + z_{P}} \right)^{2} + c^{2}}} \right)}} & (25) \end{matrix}$

Knowing that

$c^{2} = \frac{a^{2}z_{P}^{2}}{r^{2} - a^{2}}$

each term of

$z_{w}^{\prime} = {f\left( {1 - \frac{f\left( {f + z_{P}} \right)}{\left( {f + z_{P}} \right)^{2} + \frac{a^{2}z_{P}^{2}}{r^{2} - a^{2}}}} \right)}$

can be divided by z_(P) ². Thus, if z_(P)»f and (f+z_(P))≅z_(P) an approximation can be derived:

$\begin{matrix} {= {f\left( {1 - {\frac{f}{z_{P}}\left( {1 - \frac{a^{2}}{r^{2}}} \right)}} \right)}} & (26) \end{matrix}$

Since a′ is the radius of the wait of the imaged pixel beam 72, it is also the minimum value of the positive hyperbola for the imaged pixel beam 72 i.e. it corresponds to the value of x′ as defined by equation 20 for z′=z′_(W):

a′=x′(z′ _(W))  (27)

This gives:

$\begin{matrix} {{\frac{f^{2}c^{2}}{a^{2}}a^{\prime^{2}}} = {{c^{2}\left( {f - z_{w}^{\prime}} \right)}^{2} + \left( {{\left( {f + z_{P}} \right) \cdot z_{w}^{\prime}} - {f \cdot z_{P}}} \right)^{2}}} & (28) \\ {\left. \Leftrightarrow a^{\prime} \right. = {\frac{a}{f}\sqrt{\left( {f - z_{w}^{\prime}} \right)^{2} + \frac{\left( {{\left( {f + z_{P}} \right) \cdot z_{w}^{\prime}} - {f \cdot z_{P}}} \right)^{2}}{c^{2}}}}} & (29) \\ {\left. \Leftrightarrow a^{\prime} \right. = {\frac{a}{f}\sqrt{\left( {f - z_{w}^{\prime}} \right)^{2} + {\frac{r^{2} - a^{2}}{a^{2}}\left( {{\frac{f + z_{P}}{z_{P}} \cdot z_{w}^{\prime}} - f} \right)^{2}}}}} & (30) \end{matrix}$

Let's consider a point x, y and its conjugate x′, y′ through the optical system 71. For a ray of light passing by point x, y and the optical system 71 centre, the laws of optics give:

$\begin{matrix} {\frac{x}{z} = {\frac{x^{\prime}}{z^{\prime}}.}} & (31) \end{matrix}$

If the ray of light is a ray of incidence θ_(x) when hitting the optical system 71 at height x₀ then:

$\begin{matrix} {{\tan \mspace{14mu} \theta_{x}} = {\frac{x - x_{0}}{z} = {\frac{x}{z} - \frac{x_{0}}{z}}}} & (32) \end{matrix}$

and similarly:

$\begin{matrix} {{\tan \mspace{14mu} \theta_{x}^{\prime}} = {\frac{x^{\prime}}{z^{\prime}} - \frac{x_{0}}{z^{\prime}}}} & (33) \end{matrix}$

Thus:

$\begin{matrix} {{{\tan \mspace{14mu} \theta_{x}^{\prime}} = {{{\tan \mspace{14mu} \theta_{x}} + \frac{x_{0}}{z} - \frac{x_{0}}{z^{\prime}}} = {{\tan \mspace{14mu} \theta_{x}} + {x_{0}\left( {\frac{1}{z} - \frac{1}{z^{\prime}}} \right)}}}}\;} & (34) \end{matrix}$

As known from the lens maker formula:

$\frac{1}{z^{\prime}} = {\frac{1}{z} + \frac{1}{f}}$

we obtain:

$\begin{matrix} {{\tan \mspace{14mu} \theta_{x}^{\prime}} = {{\tan \mspace{14mu} \theta_{x}} - \frac{x_{0}}{f}}} & (35) \end{matrix}$

And similarly:

$\begin{matrix} {{\tan \mspace{14mu} \theta_{y}^{\prime}} = {{\tan \mspace{14mu} \theta_{y}} - \frac{y_{0}}{f}}} & (36) \end{matrix}$

Consequently, the equation representing the imaged pixel beam 72 is:

$\begin{matrix} \left\{ \begin{matrix} {x = {x_{o} + {a^{\prime}{\sqrt{1 + \frac{\left( {z - z_{w}^{\prime}} \right)^{2}}{c^{2}}} \cdot {\cos (v)}}} + {\left( {z - z_{0}} \right) \cdot {\tan \left( \theta_{x}^{\prime} \right)}}}} \\ {y = {y_{o} + {a^{\prime}{\sqrt{1 + \frac{\left( {z - z_{w}^{\prime}} \right)^{2}}{c^{2}}} \cdot {\sin (v)}}} + {\left( {z - z_{0}} \right) \cdot {\tan \left( \theta_{y}^{\prime} \right)}}}} \end{matrix} \right. & (37) \end{matrix}$

with parameters: z′_(w), θ′_(x), θ′_(y), a′ a obtained from the parameters z_(P), θ_(x), θ_(y), a of the pixel beam 70 by:

$\begin{matrix} \left\{ \begin{matrix} {z_{w}^{\prime} = {f\left( {1 - \frac{f\left( {f + z_{P}} \right)}{\left( {f + z_{P}} \right)^{2} + c^{2}}} \right)}} \\ {a^{\prime} = {\frac{a}{f}\sqrt{\left( {f - z_{w}^{\prime}} \right)^{2} + \frac{\left( {{\left( {f + z_{P}} \right) \cdot z_{w}^{\prime}} - {f \cdot z_{P}}} \right)^{2}}{c^{2}}}}} \\ {{\tan \mspace{14mu} \theta_{x}^{\prime}} = {{\tan \mspace{14mu} \theta_{x}} - \frac{x_{0}}{f}}} \\ {{\tan \mspace{14mu} \theta_{y}^{\prime}} = {{\tan \mspace{14mu} \theta_{y}} - \frac{y_{0}}{f}}} \end{matrix} \right. & (38) \end{matrix}$

with

$c^{2} = \frac{a^{2}z_{P}^{2}}{r^{2} - a^{2}}$

In a step 6010, the processor 501 of the apparatus 500 computes an intersection, in the object space of the first optical system, of a pixel beam 70 of a collection of pixel beams and a given plane, defined by its affine equation in a reference coordinate system.

Then in a step 6020, the processor 501 computes the image conjugate of this intersection through the optical system 71.

In a step 6030, the processor 501 computes the imaged pixel beam 72 from the image conjugate of the intersection computed during step 6020.

In a step 6040, the processor 501 of the apparatus 500 then computes an intersection of the imaged pixel beam 72 with the sensor of the optical system 71.

Then, in a step 6050, the processor 501 computes the integration of the light irradiance over the said intersection.

Such a method for generating data representative of pixel beams offers a precise knowledge and parametrization of a collection of pixel beams in a given focal volume. The parameters representative of the imaged pixel beams associated with RGB values associated to the pixels corresponding to said pixel beams form a parametric sampling of the focal volume of the second optical system 72 useful for image processing operations since the parametric sampling of the focal volume is homogeneous.

In an embodiment of the invention, the first optical system, which object space has been sampled by the collection of pixel beams 70 is the optical system of an actual camera whereas the second optical system 71 is the optical system of a virtual camera.

In another embodiment of the invention, the first optical system, which object space has been sampled by the collection of pixel beams 70 is the optical system of a virtual camera whereas the second optical system 71 is the optical system of an actual camera.

Although the present invention has been described hereinabove with reference to specific embodiments, the present invention is not limited to the specific embodiments, and modifications will be apparent to a skilled person in the art which lie within the scope of the present invention.

Many further modifications and variations will suggest themselves to those versed in the art upon making reference to the foregoing illustrative embodiments, which are given by way of example only and which are not intended to limit the scope of the invention, that being determined solely by the appended claims. In particular the different features from different embodiments may be interchanged, where appropriate. 

1. A computer implemented method comprising generating data representative of a volume occupied by a set of rays of light passing through a pupil of an optical system, said volume occupied by said set of rays of light being called a pixel beam; acquiring data representative of said pixel beam in an object space of a first optical system, said pixel beam being occupied by a set of rays of light passing through a pupil of said first optical system and an image conjugate, in said object space of the first optical system, of at least one pixel of a sensor associated to said first optical system; and generating data representative of said pixel beam by computing an image conjugate of said pixel beam through a second optical system.
 2. The method according to claim 1 further comprising: obtaining an intersection, in the object space of the first optical system, of said pixel beam with a plane, obtaining an image conjugate of said intersection through said second optical system to provide data representative of the pixel beam.
 3. The method according to claim 1, wherein the first optical system is an actual camera and the second optical system is a virtual camera.
 4. The method according to claim 1, wherein the first optical system is a virtual camera and the second optical system is an actual camera.
 5. The method according to claim 2, wherein the first optical system is an actual camera and the second optical system is a virtual camera.
 6. The method according to claim 2, wherein the first optical system is a virtual camera and the second optical system is an actual camera.
 7. A device comprising a light filed data generation module configured to generate data representative of a volume occupied by a set of rays of light passing through a pupil of an optical system, said volume occupied by said set of rays of light being called a pixel beam; acquire data representative of said pixel beam in an object space of a first optical system, said pixel beam being occupied by a set of rays of light passing through a pupil of said first optical system and an image conjugate, in said object space of the first optical system, of at least one pixel of a sensor associated to said first optical system; and generate data representative of said pixel beam by computing an image conjugate of said pixel beam through a second optical system.
 8. The device according to claim 7, wherein the light field data generation module is further configured to: obtain an intersection, in the object space of the first optical system, of said pixel beam with a plane, obtain an image conjugate of said intersection through said second optical system 15 to provide data representative of the pixel beam.
 9. The device according to claim 7, wherein the first optical system is an actual camera and the second optical system is a virtual camera.
 10. The device according to claim 7, wherein the first optical system is a virtual camera and the second optical system is an actual camera.
 11. The device according to claim 8, wherein the first optical system is an actual camera and the second optical system is a virtual camera.
 12. The device according to claim 8, wherein the first optical system is a virtual camera and the second optical system is an actual camera.
 13. The device of claim 8, wherein said device is configured for rendering an image from said light field data.
 14. The device of claim 7, further comprising: an array of micro lenses arranged in a regular lattice structure; a photosensor configured to capture light projected on the photosensor from the array of micro lenses, the photosensor comprising sets of pixels, each set of pixels being optically associated with a respective micro lens of the array of micro lenses; and a component for generating metadata.
 15. The device of claim 8, further comprising: an array of micro lenses arranged in a regular lattice structure; a photosensor configured to capture light projected on the photosensor from the array of micro lenses, the photosensor comprising sets of pixels, each set of pixels being optically associated with a respective micro lens of the array of micro lenses; and a component for generating metadata.
 16. A data package for data representative of a volume in an object space of a first optical system occupied by a set of rays of light passing through a pupil of said first optical system and a conjugate, in said object space of the first optical system, of at least one pixel of a sensor associated to said first optical system, said volume occupied by said set of rays of light being called a pixel beam said data representative of a pixel beam being an image conjugate of said pixel beam through a second optical system.
 17. A computer program product for a programmable apparatus, the computer program product comprising acquiring data representative of said pixel beam in an object space of a first optical system, said pixel beam being occupied by a set of rays of light passing through a pupil of said first optical system and an image conjugate, in said object space of the first optical system, of at least one pixel of a sensor associated to said first optical system; generating data representative of said pixel beam by computing an image conjugate of said pixel beam through a second optical system.
 18. The computer program product according to claim 17 further comprising: obtaining an intersection, in the object space of the first optical system, of said pixel beam with a plane, obtaining an image conjugate of said intersection through said second optical system to provide data representative of the pixel beam.
 19. The computer program product according to claim 17, wherein the first optical system is an actual camera and the second optical system is a virtual camera.
 20. The computer program product according to claim 17, wherein the first optical system is a virtual camera and the second optical system is an actual camera. 